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The first three electronically excited states in the alkaline-earth-metal atoms magnesium, calcium, 
and strontium comprise the (nsnp) ^P] (J — 0, 1, 2) fine-structure manifold. All three states are 
metastable and are of interest for optical atomic clocks as well as for cold-collision physics. An 
efficient technique — based on a physically motivated potential that models the presence of the ionic 
core — is employed to solve the Schrodinger equation for the two-electron valence shell. In this way, 
radiative lifetimes, laser-induced clock shifts, and long-range interaction parameters are calculated 
for metastable Mg, Ca, and Sr. 

PACS numbers: 31.25.Jf, 32.70.Cs, 34.20.Cf 



I. INTRODUCTION 

Effectively, alkaline-earth-metal atoms are two-electron atoms, at least at low excitation energies. Alkaline-earth- 
metal atoms display — like true two-electron systems — spin singlet and triplet level manifolds, which give rise to 
narrow intercombination transitions. However, since all alkaline-earth-metal atoms possess a closed K shell, their 
level structure does also show notable differences to that of hcliumlike systems. 

In this work we focus on the lowest "^P" manifold, with fine-structure components J = 0,1,2, in the alkaline- 
earth-metal atoms magnesium, calcium, and strontium. These three fine-structure states represent the first three 
electronically excited states in these species. The excitation energies, relative to the ^5*0 ground state, are reproduced 
in Tabled IJ. The table illustrates that for all three atoms, the excitation energies correspond to optical wavelengths 
accessible by modern laser sources. However, that does not imply the fine-structure components of the "^P" manifold 
are easy to excite. The transition from the ground state to the '^Pf is spin-forbidden, that to the ^Pq and the ^P^" 
states is, additionally, in general electric-dipole-forbidden. (See Ref. for a thorough introduction to the topic of 
forbidden transitions in one- and two-electron atoms.) This makes the '^P° manifold in the alkalinc-carth-mctal atoms 
interesting for several applications in the physics of cold atoms. 

The unique level structure triggered advances in laser-cooling technology and allowed to explore processes in a 
magneto-optical trap at a new level 0, 0, IS S S S IS ^] ■ Current successful approaches to an optical frequency 
standard u sing neutral, laser-cooled atoms are based on the relatively small linewidth of the ^5*0 ^ '^Pi° transition 
in calcium pjl 11^ . Significant improvement in accuracy is expected from the strategy of storing fermionic *^Sr 
atoms in an optical lattice and probing the ^5*0 "^Pp transition [liL lla lla llTl . 

Due to the prevalence of stable alkaline-earth-metal isotopes with vanishing nuclear angular momentum, it is possible 
to investi gate cold collisions and to perform photoassociation spectroscopy free of any comphcations from hyperfine 
structure pjl l20l l2lL 122. l23l | . However, it is not possible to form a quantum-degenerate gas |24Ll25ll2a | of ground-state 
alkaline-earth-metal atoms by employing evaporative cooling in a pure magnetic trap. Therefore, it has been suggested 
to make use of the long-lived ^P° state [27, 28]. Several groups have already succeeded in magnetically trapping ^Pf 
alkaline-earth- metal atoms [23, ISfl Is^, Is^l ■ Theoretical studies have uncovered that anisotropic interactions between 
•^Pf atoms in an external magnetic field lead to unexpectedly rich collision physics jsS 0,113 • 

It is the purpose of this paper to provide a concise overview of the properties of Mg, Ca, and Sr in the '^Pq , 
3pf , and 3P2° states. In Sec. HII we describe our specific approach to the electronic-structure problem of effective 
two-electron atoms. The method adopted is not only numerically efficient, it can supply atomic data at a level of 
accuracy comparable with high-level ab initio calculations. Section IIIII is devoted to a discussion of spontaneous 
decay of the individual fine-structure components of the ^P° manifold. Then, in Sec. IIVL we calculate the effect of 



TABLE I: The excitation energies, in eV, of the first three excited states in the alkaline-earth-metal atoms magnesium (3s3p), 
calcium (4s4p), and strontium (5s5p). The data are taken from Ref. 1]. 
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the dynamic Stark shift on the ^Sq ^ ^Pg clock transition. Parameters characterizing the long-range interaction 
between metastable alkaline-earth-metal atoms are the subject of Sec. [3 We conclude with Sec. I VII Unless indicated 
otherwise, atomic units are used throughout. 



II. TREATMENT OF THE ELECTRONIC-STRUCTURE PROBLEM 



Our strate gy t oward solving the electronic-structure problem of alkaline-earth-metal atoms involves the following 
steps [3^ • We first treat the effective one-electron eigenvalue problem of the singly-charged alkaline-earth- 

metal ion with a single valence electron. This valence electron moves in a one-particle model potential that reproduces 
the valence excitation energies of the monocation. We represent the radial degree of freedom of the electronic wave- 
function in a finite-element basis set. From the solutions of the one-electron problem, two-electron basis functions are 
constructed. The effective two-electron Hamiltonian, which describes the neutral atom and which fully incorporates 
valence-electron correlation, is represented in this basis and diagonalized. In this way, eigenenergies and eigenvectors 
of the two-electron valence shell are obtained. 



A. One-particle Hamiltonian 



When treating the effective one-electron problem of the monocation, we employ the one-electron Hamiltonian 

hi ^ -^V^ + V + V'^'°\ (1) 

The operator V represents the electrostatic field generated by the noble-gaslike core. This is a central field. However, 
an explicit dependence of the associated potential on the orbital angular momentum quantum number I must be 
included to take account of the fact that electrons with I — dive deep into the core and experience relativistic effects 
(orbital contraction). Such scalar relativistic effects decrease with the average distance from the nucleus and thus 
with I. Therefore, we assume that the spin-angular representation of V can be written as 

{j,mj,l,s\V\j',m.i>,l',s) = (5jj'(5„^,m^.,(5i,i'V/(r). (2) 

The quantum numbers j and rrij refer to the total angular momentum of the valence electron, s = 1/2 denotes its 
spin. In this work, the following parametrization of Vi{r) is used: 



Vi{r) = {2 + (Z - 2) exp (-a/,ir) -I- a;,2''exp {-ai^sr)} 
r 



2r4 ^ 



exp I 



(3) 



Z is the nuclear charge. The parameters ai^i, ai^2^ o/,3i c«cpj and r; are listed in Ref. |3q . The potential V;(r) describes 
the interaction of the valence electron with the ionic core at various length scales. For very large distances from the 
core, the electron is attracted by a point charge of charge -1-2. As the electron comes closer, the ionic core responds to 
the presence of the electron and becomes polarized, as expressed by the term proportional to a^plr'^. Below r = ri, 
the electron dives into the core. The parameters a;,i, a/, 2, and a/, 3 mediate the transition from the exterior region of 
the core to the interior, where at very small length scales the electron interacts with the unscreened charge Z. 
Spin-orbit interaction is represented in Eq. by the operator V^^°^: 



where 



%^f{r) = {j{3 + l)-l{l + l)-s{s + l)} 



(4) 
(5) 
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and a w 1/137.036 is the fine-structure constant. 
The eigenstates of hi can now be written as 



\j,mj,l,s) 



(6) 
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and we can focus on solving the radial equation 



(7) 



To this end, we apply a technique based on finite elements js^ liolllll lil. li^lli^l . 



(8) 



B. Finite-element basis 

For the radial degree of freedom, we introduce a quadratically spaced grid of A'^ + 1 grid points between r = and 
= Tinax, i-e., Ti = ^max*^/-^^, * = 0, ...,A^. In each interval [ri,ri+i], six unique, linearly independent fifth-order 
Hermite interpolating polynomials can be constructed that satisfy the boundary conditions 



dr'^ 



dr^ 



0, 
0, 



(9) 



where j, fc = 0, 1, 2. The functions are defined to vanish everywhere outside [r^, ^i+i]. Using the definitions 

3o(r) = /(r)3 [6/(r)2 - 15/(r) + 10] , 
31 (r) = Ar/(r)3[3/(r)2-7/(r)+4], 



(10) 



we have, for r e [r^, r^+i] and j = 0, 1, 2, 



and 



Ar 



Ar 



(11) 



(12) 



These formulas correct some misprints in Eq. (49) of Rescigno et al. 43] . Three continuous basis functions, with 
continuous first and second derivatives, can now be associated with grid point (i = 1, N — 1): 



BjAr) = Pj.(r) +Q,,,-i(r), j = 0,1,2. 



(13) 



Using these, the function value at grid point r^, the first derivative, and the second derivative of any wavefunction 
can be represented. At the end points, and r^r, we set 



Bo,o(r)-0, Bj,o{r) = P,,o{r), J = l,2, 
which is consistent with the boundary condition specified in Eq. and 

BaM{r)^0, S,- Ar(r) = g,, jv-i(r), j = l,2, 
which selects solutions to Eq. ((TJ that vanish at r = Tmax- 



(14) 



(15) 
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The solutions to the eigenvalue problem of h, 

huk{r) = EkUkir), (16) 

(the explicit dependence on angular momentum quantum numbers is suppressed) are found by expanding the eigen- 
fimctions in terms of the finite-element basis functions Bj^i{r): 

N 2 

Hence, we have to solve the generalized eigenvalue problem 

hf3k = ekOf3k, (18) 

where 

(^)o,0.O',<') ^ / '"^'^ drB,4r)hB,,,Ar), (19) 



= I drB,4r)By^e[r), (20) 

and 

(^fe)(j,i) = (21) 

Since the basis function Bj i{r) vanishes outside [r^-i, ri+i], both matrices, h and o, have a simple banded structure 
with small bandwidth. The high degree of sparsity is ideal for iterative solvers ■ Using the Lanczos-based package 
ARPACK |46j . we calculate — for each valid combination of angular momentum quantum numbers I and j — the first 
18 eigenfunctions outside the core shells. In other words, solutions associated with inner shells are skipped. Thus, 
the selected valence-electron solutions display correct nodal behavior. 



C. Two-particle Hamiltonian 



We describe the neutral atom employing the effective two-electron Hamiltonian 

h2 ^ hi{l) + hi{2) + i (22) 

\Xi - X2\ 

i.e., within the approach we take, both valence electrons move in the field of the model potential V (Eqs. ^ and lO), 
are subject to spin-orbit coupling through 'V^^°^ (Eqs. Q and lO), and experience mutual Coulomb repulsion. 

In order to tackle the eigenvalue problem of /i2, we construct two-electron basis functions, in H coupling, from the 
solutions of the one-electron problem: 



$ = -4^^^1^i^^^^^^^ E C^(jiJ2J;™,,m,,Af)b-i,m,,,;i,si)|j2,m,,,;2,S2), (23) 

where A is an antisymmetrization operator and C{jij2J]Tnj-^mj^M) is a Clebsch-Gordan coefficient. The total 
electronic angular momentum, J, the associated projection, M, as well as the atomic parity, H = (— 1)'!+'^, are the 
only good quantum numbers in the two-electron problem, Eq. H22(l . The matrix representation of /ii(l) -f hi(2) — the 
Hamiltonian of two noninteracting valence electrons — in the basis of the two-electron functions is, of course, diagonal. 
Using the familiar relation jl^l 



1 \ - X - Ait 



^jTT^^yiM.Vi)yM.V2), (24) 



l—O m— — l ^ 
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r< = min{ri,r2}, r> = max{ri,r2}, 

together with standard Wigner-Racah algebra [isL l49l fsof , the matrix of the electron-electron interaction term in 
Eq. 122|) can be calculated. 

The resulting matrix representation of /i2 is sufficiently small to allow for full numerical diagonalization. Properties 
of the valence shell can now be computed from the spectrum and eigenvectors of this matrix. In order to estimate the 
accuracy of our calculations, we first calculate the physical quantity of interest (a radiative decay rate, for example), 
ensuring convergence with respect to all basis set parameters. In this step, the eigenenergies and the relevant matrix 
elements with respect to the eigenvectors are taken from our calculations. Then we recalculate that physical quantity 
by replacing the energies of the most important states by experimental energies taken from Ref . jlj] . Finally, we quote 
the result based on the experimental energies together with an error estimate given by the difference from the result 
calculated using theoretical energies. 



III. RADIATIVE DECAY RATES 



3 po 
M3 



The bosonic atoms ^"^Mg, '^'^Ca, and ^^Sr are by far the most abundant isotopes of magnesium, calcium, and 
strontium, respectively. The nuclei in these isotopes have vanishing nuclear spin, implying that the quantum number 
J = in the state is exact. As all one-photon operators have a rank greater than zero, the most important 
way the ^Pq state can decay is via an ElMl two-photon process, which, in contrast to ElEl, can mediate the parity 
change between ^Pq and \Sn- The ElMl process has been analyzed in detail for the radiative decay of the ^Pq state 
in heavy heliumlike ions |53,|52|, where this process competes with the usual El decay. (In heliumlike systems, the 
^Pq state is not the first excited state.) 

The ElMl-decay rate from the ^Pq state, symbolized by the state vector \i) with energy Ei^ to the ^5*0 ground 
state (I/), Ef) can be readily derived within the framework of quantum electrodynamics !53| and is given by 



ElMl 



277r 



"'0 



dL02^2 



E 



(/ II Ml II n+) (n+ II D \\ i) 



E 



En. + UJi — Ei 



(/ II D II n_) (n_ II Ml II z) 



En_ +UJ2- E, 



S{Ef +LUi+LU2~ E,). 



(25) 



We use the notation (• || • || •) to denote a reduced matrix element |4^^, "S^l. Equation (|25|l illustrates that the two 
emitted photons with energies loi and UJ2, respectively, share the available energy, Ei — Ef . The decay can proceed in 
two indistinguishable ways. In one of the two, the El photon is emitted first, mediated by the electric-dipole operator 



D, = -Y,nCi,g{'d,,^^), (26) 

i 

and a virtual transition to the manifold of states \n^) with even parity and J — 1 takes place. (We use the definition 
C;,q = ^/4n/2l + lYi^q.) The intermediate states are then coupled to the ground state by the magnetic-dipole operator 



Ml.q = -flB (J, + Sq) , (27) 

fj,B denoting the Bohr magneton. Alternatively, we first have a virtual Ml transition to states |rt_) with odd parity 
and J = 1, followed by El emission. In the case of the lowest "^Pp state in Mg, Ca, and Sr, the energy denominators 
in Eq. (|25|l differ from zero for all uji and 1^2, so that there are no numerical complications arising from poles. 

Our calculated ElMl rates are presented in Tabled They are, not surprisingly, very small throughout. Even 
in strontium, which decays faster by more than an order of magnitude than calcium and magnesium, the lifetime is 
about 5800 years. This lifetime is comparable to the half-life of the radioactive isotope "'^^C. 
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TABLE II: Spontaneous decay rate, in s^^, of the '^Pg state in the alkahne-earth-metal atoms magnesium (3s3p), calcium 
(4s4p), and strontium (5s5p). The spinless bosonic isotopes can decay only via ElMl two-photon emission. The stable 
fermionic isotopes — ^^Mg, *^Ca, and *^Sr — decay by a hyperfine-induced El one-photon process. 

Mg Ca Sr 

boson 1.6(1) X 10"^^ 3.9(2) x 10"^^ 5.5(3) x 10"^^ 
fermion 9(3) x 10"^ 3(1) x 10''^ 9(3) x 10'^ 



The nticlear angular momentum, 7, in the fermionic isotopes of magnesium, calcium, and strontium differs from 
zero [13: /(^^Mg) = 5/2, iC^^Ga) = 7/2, /(**'^Sr) = 9/2. The finite nuclear magnetic-dipole moment couples the ^Pg 
state to states of the same parity but with J = 1. The most important among these states are the energetically lowest 
^Pf and "'^Pi states. These can decay directly to the ground state by emission of an El photon. 

The interaction of the valence electrons with the magnetic-dipole moment, M]™'^ , and the electric-quadrupole 
moment, Q™*^, of the atomic nucleus is given by the hyperfine-coupling operator |55l | 



We employ the notation 



B 



k2\K,M 



E 



C{kik2K; mim2M)Ak^,m^Bk2^rn2 



(28) 



(29) 



for the rank-X tensor product of tensors (rank ki) and B^^ (rank k2). 

It is not a priori clear that the valence-electron wavefunctions derived from the model potential, Eq. (Q, have the 
correct behavior near the nucleus, even if valence excitation spectra can be well reproduced. In order to check this, we 
employed the hyperfine-coupling operator Vm to quantitatively calculate hyperfine parameters for the (5s5p) and 
^P2 states in ^^Sr 54). For ^Pf we find A = -278 MHz and B = -30.8 MHz. This result is in reasonable agreement 
with an experiment by zu Putlitz |56[, who measured A — -260 MHz and B — -35.7 MHz. Heider and Brink [s^ 
determined the hyperfine parameters for the state: A = -213 MHz and B — 67.2 MHz. Our calculation gives A 
= -231 MHz and B = 56.5 MHz. Even though the agreement is not perfect, it indicates that we should be able to 
calculate the hyperfine-induced coupling of the ^P§ state to other electronic levels with an accuracy of about 10%. 

Since the hyperfine-coupling operator Vm (Eq. H28|l ) is a tensor of rank with respect to rotations of nucleus and 
electrons, it conserves the total angular momentum of the atom, F. Furthermore, for the transition from ^Pq to ^5*0, 
F=I. First-order perturbation theory leads to the following formula for the hyperfine-induced spontaneous decay rate: 



Fhf = -- 



4 a^uj^ 



9 2/ -hi 



(30) 



^ (/ \\D\\ n-){F,n 



Vm II F,i) 



En 



Here, uj = Ei — E f is the energy of the emitted photon, and as before, the states 



have odd parity and J = 1. 

The electric-quadrupole term in Vm does not play a role, for it couples "^Pq to states with J = 2. These, however, 
cannot decay to the ground state via El one-photon emission. We use the abbreviated notation |P, n) in the reduced 
matrix elements of Vm in Eq. H30(l to indicate that these states are angular-momentum coupled states of the entire 
atom with total angular momentum F. 

The ■^Pg decay rates in the fermionic alkaline-earth-metal isotopes, as calculated on the basis of Eq. (|30l) . are 
collected in Table We note that the lifetimes are now much shorter — by about 10 orders of magnitude — than 
in the bosonic isotopes. Nevertheless, this is still sufficiently long to make the ^Sq ^Pq transition interesting for 
optical atomic clocks. Especially the strontium isotope *^Sr — which has a rather narrow ^Pq linewidth of about 1 
mHz — is a candidate for a future ultraprecise atomic clock In fact, this is the only one of the three isotopes 

25Mg, 43Ca, and **^Sr, for which we are aware of an analogous '^Pq decay rate determination: Pal'chikov |58j | obtained 
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TABLE III: Spontaneous decay rate, in s ^, of the state in the alkahne-earth-metal atoms magnesium (3s3p), calcium 
(4s4p), and strontium (5s5p). The decay is dominated by the El rate for the transition ^P" ^So. 

Mg Ca Sr 

3.6(1) X 10^ 2.1(2) X 10^ 4.1(4) x 10"^ 



a rate of 5.5 x lO^'^ s"^ for this ®^Sr transition. This is somewhat smaUer than our result, apparently for two reasons. 
First, we have used experimental excitation energies, and second, we have included a larger number of intermediate 
states when evaluating Eq. I|30|l . If we apply the same restrictions as Pal'chikov, we find 6 x 10^"^ s^^. 



B. =^Pi° 

The '^P° state can decay via two obvious decay modes. On one hand. Ml one-photon decay from the '^P" to the 
•^Pq state can take place. This mode, however, is slower by many orders of magnitude than the El decay directly to 
the ^5*0 ground state. The only reason the state is metastable is the necessity of a spin flip. The El decay is 
enabled by the presence of spin-orbit coupling — a relativistic effect that becomes more pronounced in heavier atoms. 

The El decay rate is given by 0| 

rE.^3^|(/||z.|M)r (31) 

Our calculated reduced electric-dipole matrix elements for the '^P° ^ ^Sq transition agree to about the 10% level with 
recent high-level many-body and configuration-interaction calculations 59, 60] . The result for the calculated decay 
rate, Fei, is presented in Table ITTTl As expected, the decay rate increases from Mg to Sr due to enhanced spin-orbit 
coupling. 

The '^Pi decay rate in all three alkaline-earth-mctal atoms is sufficiently large to allow their direct experimental 
determination. The most recent measurements give 1.9(3) x 10^ s~^ for Mg [bj, 2.9(2) x 10*^ s~^ for Ca [b^, and 
4.7(1) X 10* s~^ for Sr [b^. Perfect agreement is not expected for our simple model approach. Nonetheless, we may 
conclude that our calculated reduced electric-dipole matrix elements are apparently accurate to about the 10% level. 



Alkaline-earth-metal atoms in the "^Pq state do not possess a magnetic-dipole moment. In the ^P° state they do, 
but even in Mg this state decays after about 5 ms. The third candidate, ^P2, is the only one useful for magnetic 
trapping [g^I, due to its inherently long lifetime p7| . 

There are several different modes through which the '^P^' state can decay. Transitions within the "^P" fine-structure 
manifold are mediated by parity-conserving one-photon operators. The rate associated with a magnetic-dipole tran- 
sition to the ^P° state can be written as |53 

rMi = ^^^|(/||Mi||z)r (32) 

While it is permissible to neglect magnetic-octupole decay, the electric-quadrupole decay to the ■^Pf state must be 
taken into consideration. The rate for the latter process is given by 55] 

where 

Qg = -J2r^C2^M,V^^)■ (34) 

i 

The E2 mechanism also allows for decay from ■^Pa" to ^Pq. Finally, direct decay to the ^Snground state requires a 
parity-changing one-photon operator of rank 2, i.e., the magnetic-quadrupole operator jB^. l65l] 
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TABLE IV: Spontaneous decay rates, in s ^, of the ''P2 state in the alkaline-earth-metal atoms magnesium (3s3p), calcium 
(4s4p), and strontium (5s5p). 



+-rQ n CI i" ion 


type 


Mo- 


Ca 


Sr 




Ml 


1(1) X lO"'' 


1.6(3) X 10"^ 


8(2) X 10-* 


^ p° ^ 


E2 


1(3) X 10"^^ 


3(1) X 10-1" 


3(1) X 10"^ 


3 po ^ 3 po 

-^2 ^ -M3 


E2 


0(1) X 10'" 


9(3) X 10-1" 


9(4) X 10"^ 




M2 


4.2(2) X 10""* 


1.1(1) X 10-'' 


1.1(2) X 10"* 


total rate 




4.2(2) X lO'"* 


1.3(1) X 10-* 


9(2) X 10-* 


Ref. [27] 




4.42 X 10"" 


1.41 X 10"" 


9.55 X 10-* 



M2,q -V&a r, I [Ci(i9„ ip.,) ® s,]2,g + ^ [Ci(z9^, (/?,) ® h]2,q \ ■ (35) 
The associated rate is ^91 



M2 



15 2 J, + 1 



l(/l|A^2|M)|'. (36) 



Our calculated decay rates are shown in Table HVl We see that in Mg and Ca, the exotic M2 mode dominates by 
far. In Sr, however, due to the uJ^ dependence in Eq. 1)3211 . Ml decay makes the biggest contribution. Also shown in 
Table HVl is a comparison of our calculated total rates with those determined by Derevianko, Ref. We find good 
agreement and conclude that the ^Pj" hfetimes in Mg, Ca, and Sr are, respectively, 2.4(1) x lO'^ s, 7.7(6) x 10^ s, and 
1.1(3) X 10'^ s. The first measurement, by Yasuda and Katori, of the lifetime in ^^Sr yielded 5001^3q s j6^, which 
is compatible with the current level of theory. 



IV. DYNAMIC STARK EFFECT 



As mentioned in the introduction. Sec. HI the transition from the ground state to the ^Pq state in the fermionic 
isotope ^^Sr is of interest for the development of an ultraprecise optical clock based on neutral atoms trapped in an 
optical lattice j3 lisl EH IitI- 0| • It has been pointed out by Katori et al. that with this approach a relative accuracy 
of better than 10"^'' should be achievable [l^ . 

The optical lattice serves three purposes. If each lattice site is occupied by at most one atom, collisional clock 
shifts can be suppressed (see Ref. [l^ for a discussion of frequency shifts caused by dipole-dipole interactions among 
atoms in different lattice sites). Furthermore, using trapped atoms, long interaction times with the probe laser can be 
realized. This is particularly important when exciting a transition with exception ally small linewidth (see Table ITHl . 
A third benefit is the possibility of confining the atoms in the Lamb-Dicke regime [ij, [1^, [l3i US ■ In this way, atomic 
recoil- and Doppler-shifts can be eliminated. 

The intense laser field needed to confine the atoms perturbs them. To lowest order in perturbation theory, each 
atomic level is shifted by a level-specific amount proportional to the intensity of the laser |63|- This effect is referred 
to as dynamic (or ac) Stark shift. The dynamic Stark shift of state |m) can be calculated from [68l | 

. ^ 27r aX ,„_, 
= — , . (37) 

\{n\\D\\m)\' ^ \{n\\D\\ 

Em — En — LO E,n — En + LO \ 

In this equation, I is the laser intensity in units of Xq = e^/ag = 6.4364 x 10^^ W/cm^; w is the laser photon energy. 

However, it is possible to find a magic wavelength for the t rapping laser field such that the dynamic Stark shifts 
of the two levels involved in the clock transition are identical [l4. Il6l IT^ . In other words, at the magic wavelength, 
the clock frequency is virtually unaffected by the presence of the optical lattice. In a recent calculation, the magic 
wavelength for the ^S n — » ^Pq transition in Sr was found to be about 800 nm T^, in agreement with the experimental 
value of 813.5(9) nm |l3. 
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TABLE V: Magic wavelength, in nm, for the transition (ns^) ^So (nsnp) ^Pq. At this wavelength, the dynamic Stark shifts, 
Eq. H37|l . of both levels, ^So and ^Pq, are identical. 

Mg Ca Sr 
470(f0) 700(50) 770(60) 



TABLE VI: Ce long-range dispersion coefficient, in atomic units, for two (ns^) ^So alkaline-earth-metal atoms. 

Mg Ca Sr 

This work 620(5) 2f50(60) 3260(f00) 

ReL [75] 627(f2) 222f(f5) 3170(196) 

Ref. ^ 648 2042 3212 



Using Eq. (j^H), we determined the magic wavelength for the proposed clock transition within our model-potential 
approach. The results are collected in Table IVl Within our error bars, our result for Sr is in agreement with the 
previous calculation and with the experiment. To our knowledge, Mg and Ca have not been considered before. As 
can be seen in Tabled the magic wavelength in Mg is too blue to be ideal for trapping with current laser sources. 
The magic wavelength we find in Ca appears to be much more attractive. In view of this result and of the small decay 
width of the state in ^'^Ca (see Table Hljl . fermionic calcium appears to be an excellent candidate for an optical 
atomic clock and may serve as an interesting alternative to ^^Sr. 



V. LONG-RANGE PARAMETERS 
A. 



While collisions among atoms are to be minimized in atomic clocks, interatomic interactions at temperatures close 
to absolute zero are highly intriguing and are responsible for some of the most interesting effects in the physics of 
ultracold gases .fiSj 70j 71] . In order to understand cold colUsions, detailed knowledge about the long-range 
behavior of the interatomic interaction potential is indispensable |69| . 

Here, we concentrate on interactions between two identical alkaline-earth-metal atoms in the same quantum state. 
If the two interacting atoms have vanishing an gula r momentum, J = 0, then at the leading, second-order level of the 
long-range expansion of the interaction energy [7^ l73l '74'] , both atoms can only become polarized — and thus reduce 
the overall energy of the system — by making virtual electric-dipole transitions to states with J = 1. In this case, the 
electric dipole-dipole dispersion interaction energy between atoms in atomic state |m), separated by distance R, is 
well known to be given by the leading -Ce/R^ term, where [tUtHtII 

^ _ 2 ^ Km \\D\\ m)f\{n2 \\D\\ m)f 

For ^Sq atoms, the states |?^2) have angular momentum J 1 and odd parity. The Cg coefficients we calculated 
for Mg, Ca, and Sr in the atomic ground state are listed in Table I VII Our numbers can be compared with high-level 
ab initio data from Refs. [lEl and [l^- The high quality of our results is evident. The monotonic increase of the 
Cg coefficient from Mg to Sr is a simple consequence of the greater polarizability of the heavier alkaline-earth-metal 
atoms. 



More interesting for future experiments with metastable atoms is the long-range physics of the "^Pq state. As we 
have seen in Table Ull in this state the lifetime, under spontaneous emission, of bosonic alkaline-earth-metal atoms is 
virtually unlimited. Any quenching of the metastable state will be induced by collisions with either atoms or photons. 
It is not clear, for example, on what timescale a gas consisting of '^Pq atoms will disintegrate due to inelastic processes, 
how these processes can be controlled, or what role elastic collisions have to play. An important first step toward an 
answer consists again in determining the corresponding Cg, coefficient. For -^Pq atoms, the states |n2) in Eq. lp!H| 



10 



TABLE VII: Cq long-range dispersion coefficient, in atomic units, for two (nsnp) ^Pq alkaline-earth-metal atoms. 

Mg Ca Sr 

980(30) 3020(200) 5260(500) 



TABLE VIII: Atomic electric-quadrupole moment, in atomic units, in the (nsnp) state. 

Mg Ca Sr 

This work 8.38 12.7 15.4 

Ref. [33j 8.46(8) 12.9(4) 15.6(5) 
Ref. [27] 8.59 13.6 16.4 



have angular momentum J = 1 and even parity. The result of our calculation is shown in Table IVTll As for ground 
state atoms, we observe a monotonic dependence of the Cq coefficient on the atomic number. Additionally, it should 
be noted that in the '^Pg state the Cg coefficients for all three atomic species are greater than in the ground state by 
a factor of about 1.5. This expresses the fact that for excited atoms, the relevant energy denominators in Eq. H38|) 
are smaller on average. 



As Derevianko pointed out the very long lifetime (see Sec. IIII Cll of alkaline-earth-metal atoms in the lowest 
'^P2 state makes magnetic trapping of these metastable species experimentally practical. This expectation has turned 
out to be correct |29l l3Ct l3ll l32ll . However, what ultimately limit the efficiency of evaporative cooling in a magnetic 
trap are inelastic collisions '391, not radiative decay processes. The strong coupling to inelastic scattering channels is 
caused by the anisotropic long-range interactions between the nonspherical ■^Pj" atoms. 

The most important parameter in this context is the atomic electric-quadrupole moment |27t l33 | . For an atom in 
state |m) with angular momentum J„j, the electric-quadrupole moment can be written as |34| 



Q^^J- Jrn{2J,n 1) {m \\ Q \\ m) . (39) 

Q vanishes, of course, if Jm = or Jm = 1/2. The atomic electric-quadrupole moment gives rise to direct, first- 
order quadrupole-quadrupole interactions proportional to jT^]. Derevianko was the first to calculate Q for ^Pf 
alkaline-earth- metal atoms |23; improved data were published recently We list the results from Refs. [23| and 
[33| together with our calculated values in Table IVIIII Again, our model calculations provide data of quality similar 
to the data derived from high-level many-body theory, Ref. [33| . The monotonic increase of Q with respect to atomic 
number reflects the increase in spatial extension of the excited np orbital, where n = 3 for Mg, n — 4 for Ca, and 
n = 5 for Sr. 

The electric dipole-dipole dispersion physics is also more complicated for J = 2 than for J = 0. We introduce 
intermediate dispersion coefficients to characterize the dispersion interaction between nonspherical atoms |34|: 

^ y- Kni lie II m)f\{n2 \l D \l m)f 

ni,n2 

In our specific case of two interacting '^Pj atoms, Jm — 2, and \ni), |rt2) denote even-parity states with angular 
momenta Ji and J2, respectively. Both Ji and J2 can vary between 1,2, and 3, so there are six nontrivial intermediate 
dispersion coefficients. They are displayed in Table Hxl (The remaining three coefficients can be obtained by simple 
sign changes, according to Eq. (|40|l .) Derevianko et al. use intermediate dispersion coefficients defined as 

27 

C;i''' ^-^{2J+l)^Bj,j,. (41) 



Their converted ab initio data are also shown in Table Hxl The agreement between both data sets is satisfying. 



11 



TABLE IX: Intermediate dispersion coefficients for two interacting alkaline-earth-metal atoms in the (nsnp) '^P^ state. See 
Eq. I|4()|l for a definition of Bj^ . The results presented in Ref. have been converted using Eq. JUJ. AH data are in atomic 
units. 







Mg 


Ca 


Sr 


^1,1 


(this work) 


-35.6(2) 


-81(3) 


-139(7) 


Bi,i 


(Ref [33]) 


-38(4) 


-92(9) 


-158(16) 


B2,l 


(this work) 


42.5(2) 


119(5) 


196(9) 


B2,l 


(Ref. [33]) 


44(4) 


123(12) 


203(20) 


B2,2 


(this work) 


-51.9(2) 


-176(8) 


-280(10) 


B2,2 


(Ref [33]) 


-52(5) 


-167(17) 


-264(26) 


Bs,! 


(this work) 


-73.4(7) 


-203(10) 


-370(30) 


Bs,! 


(Ref. [33]) 


-77(8) 


-225(23) 


-415(42) 


S3, 2 


(this work) 


88.6(7) 


302(20) 


546(50) 


^3,2 


(Ref [33]) 


90(9) 


306(31) 


555(56) 


Bs.s 


(this work) 


-152(2) 


-553(70) 


-1210(200) 


S3, 3 


(Ref. [33]) 


-156(16) 


-600(60) 


-1290(130) 



The restriction to long-range expansion terms up to 1/ is valid only for interatomic distances of the order of 
100 Bohr radii or larger. In this regime, however, energies associated with rotations of the diatomic frame of an 
interacting pair are comparable with the interatomic interaction energies, and the coupling between the rotational 
and the electronic d egre es of freedom induced by the anisotropic interatomic interaction must be explicitly taken into 
account. Reference |34j develops a tensorial theory appropriate for this purpose and shows how the parameters Q 
and Bj-^j^ enter the theory. 

VI. CONCLUSION 

In this paper, we have theoretically investigated properties of the energetically lowest manifold in the alkaline- 
earth-metal atoms Mg, Ca, and Sr. The latter two, in particular, are the subject of strong current interest. Using an 
electron-correlation approach for the valence shell in combination with a physically well-motivated potential for the 
valence interaction with the core shells, we have been able to reproduce existing ab initio data with good accuracy. 
This provides an independent check on a number of crucial atomic parameters for present-day experiments. 

The calculations presented here also supply data that have not, to our knowledge, been presented elsewhere. 
First, the ElMl decay rate of the ^Pq state in spinless alkaline-earth-metal species was not previously known. The 
difference in level structure to heliumlike atoms may motivate a direct search for the exotic ElMl decay process in 
heavy magnesiumlike (or berylliumlike) ions. Second, our calculation of the hyperfine-induced El decay rate of the 
•^Pq state in ''^Ca, and also our determination of the magic wavelength for the corresponding clock transition from 
the ground state of calcium, may be of interest for the development of a future generation of optical atomic clocks. 
Finally, we have calculated the Cg dispersion coefficient describing the long-range interaction between two metastable 
■^Pg alkaline-earth-metal atoms. This system, which is free from any losses due to Penning ionization, may be ideal 
for studying cold collisions between spherical, metastable atoms. 
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